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In order to assess the role of ghosts in cosmology, we study the evolution of linear cosmological 
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vector perturbations can no longer be ignored and that scalar modes diverge in the newtonian gauge 
but remain bounded in the comoving slicing. 
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I. INTRODUCTION 

Consider the action 

S=^Jd^x^R-^Jd^xV^ (5^0 d^4> + 2Vm C^,p„ C^'^p" , (1.1) 

where g is the determinant of the metric g^^, i? is the scalar curvature and Cf^^pa is the Weyl tensor. -'^ 

The two first terms describe Einstein's gravity minimally coupled to a scalar field with potential 1^(0) (and can 
also be seen as the Einstein frame formulation of /(i?) theories of gravity, see |l!]). The last term was first introduced 
by Weyl (see Q for a review of the early literature) and has ever since been present on the market of gravity 
theories, either, in recent decades, as a quantum correction popping up from various theories of quantum gravity 
(starting with [J]) or, more recently, as a phenomenological modification of Einstein's General Relativity to account 
for e.g. dark matter or energy, see e.g. [![. 

Extremisation of (|l.ip with respect to the metric yields the equations of motion: 

^ „ ^ / T^. = 9.0-5^. (ia.^a^ + n^)) , 

- Bf,y = KTf,y with i \Z J (1.2) 

where G^i, is the Einstein tensor and Bpy is the Bach tensor The divergence of the left-hand-side being identically 
zero (generalized Bianchi identity), the Klein-Gordon equation for the scalar field is redundant. The structure of this 
action and equations of motion has been thoroughly studied, in particular their hamiltonian formulation, see [tJ , as 
well as some of their solutions, for example those of Bianchi type I, see Q. 

Equations (|1.2p are fourth order differential equations for the metric components. They therefore possess extra, 
"run-away", solutions compared to the Einstein, 7 = 0, ones, which can drastically modify the predictions, even in 
the small 7 limit. 

Indeed, as shown in [9|], the theory possesses ghosts when linearised around Minkowski spacetime, that is, the 
hamiltonian contains negative kinetic terms and, as a consequence, the energy spectrum of the metric perturbations 
is not bounded from below (just as in the toy model with lagrangian (D^)^ studied earlier by Pais and Uhlenbeck, 
[iQi]). In fact most "higher derivative theories", that is, yielding equations of motion of differential order higher than 
two, are thought to possess ghosts (to the notable exception of f{R) theories of gravity, see [i|). 

Although the presence of these ghost degrees of freedom is harmless at linear level around Minkowski spacetime on 
which all modes propagate independently of each other (see e.g. jTlJ), there are strong arguments to predict that they 
yield a catastrophic collapse of any system when coupled to other fields, their energy running down to minus infinity 
in a finite time (see e.g. [lol fl^l. However, since the introduction of coupling implies that the equations of motion 
become non-linear, this catastrophic behaviour has been explicitly exhibited on toy models only, see e.g. [1311 . By the 
same token, most proposals to tame ghosts have been also illustrated by toy models only, see e.g. [l3, Il5l|. Showing 
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explicitly how the ghosts present in the particular theory of gravity described by (|l.lH1.2p may render it unviable 
when self-coupling or coupling to external fields is introduced, has not been done so far. 

Now it may happen that the malignancy of ghosts shows up already at linear level, if the background is richer than 
Minkowski spacetime. However little has been done in this direction. In ^TE\ the Hawking-Hertog proposal [14] was 
used to tame the tensor ghosts on a de Sitter background. In fvi\ the equations of motion for the scalar perturbations 
on a Friedmann-Lemaitre background were spelt out but not thoroughly analysed. 

In this paper, we aim at assessing the role of the Weyl term on the evolution of linear cosmological perturbations 
when the Friedmann-Lemaitre background is that of single-field inflation. After having obtained the equations of 
motion for the perturbations, as well as the action from which they derive, we analyse the evolution of the modes. 
We see that tensor modes are not drastically modified by the presence of ghosts. Vector modes (which have been so 
far ignored in the literature) are no longer absent as in Einstein's theory but do propagate when the Weyl term is 
present. Finally we give a master equation for the evolution of all scalar modes and find that their evolution is highly 
gauge dependent: they are unstable in the newtonian gauge but decay in the comoving slicing. We conclude on what 
should be done next. 



II. COSMOLOGICAL BACKGROUND 

When the metric is that of a conformally flat Friedmann-Lemaitre spacetime, ds^ — a{rj)^ {—drj^ + dx'^), the Weyl 
term does not contribute and the equations of motion (|1.2D for the scale factor 0(77) and the background inflaton 
(/) = f{ri) are: 

Iv''^ =H^-H', Ko^V = 2%^ +n' , (2.1) 

where a prime denotes differentiation with respect to conformal time 77 and where "H = a' /a. The Klein-Gordon 
equation which entails, 

(p" + 2-H -I- = , (2.2) 

where V^^ = ^\^=^, is also useful. 

When V{(f>) oc^" these equations are those of chaotic inflation flS\ and a detailed analysis of their solutions can 
be found in e.g. [19|. In a nutshell: after a transitory period the scale factor increases quasi-exponentially in cosmic 
time t — p a dr], while the scalar field slowly decreases. At the end of inflation the scalar field oscillates and settles 
at the bottom of its potential well while the scale factor increases in average as some power of t (t^^^ if V{(l)) cx 0^). 

III. EQUATIONS OF MOTION 

Consider the perturbed metric 

ds^ = a(77)2 [-(1 + 2A) drf + 2B, dx' drj + ((5^ + h,j) dx' dx^] (3.1) 
and perform the scalar-vector-tensor decomposition ^2^ 

B,^d,B + i3,, h,j = 2C % -f 2ay E + d^Ej + d^E, + hj (3.2) 
with diB'^ ~ diE'^ = dji^^ — Ji^ = 0. In an infinitesimal coordinate transformation the following six quantities 

■^n = A + n{B- E') + {B- E'Y , ^n = C + n{B-E'), ^,^B,-El, (3.3) 
are invariant [l^]. As for the perturbation 54> of the scalar field it is such that 

Xn = 5cl> + ^' {B-E') (3.4) 

is gauge invariant. 

We choose to work in the coordinate system B = E = 0, Ei = ("newtonian" or "longitudinal" gauge [l^, hence 

the subscript n appended to the gauge invariant scalar perturbations), so that the perturbed metric and scalar field 
reduce to 

ds^ = a{-qf {-(1 + 2*„) dif + 2*, dr^dx' + [(1 + 2$n) + hj] dx' dx^) , 5(j) = Xn- (3.5) 
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The necessary ingredients to expand the equations of motion (|1.2p at hnear order in the perturbations are given in 
Appendix [21 The resuh is: 



and 



6H $n - 2A$n - « x'n + 2a' V *„ + Xn) 



3a2 



Xn ^ _7_ 

2 ~ 3a2 



AW' . 



(3.6) 
(3.7) 



(3.8) 



where = - *n, □ = ?7^'' ^m"' and A = S'i d,j. 

These are seven equations for seven unknown quantities, two, Eq. (j3.6D . for the two components of the tensor 
perturbations hij, two, Eq. p.7p . for the two components of the vector perturbations and three, Eq. p.Sp . for the 
two scalar perturbations of the metric, and and the perturbation Xn of the scalar field. The first Eq. (|3.8p 
is the (OO)-component of the scalar part of the field equations, the second is their (Oi)-components, the third is the 
(ij)-components {i 7^ j) of their spatial part.^ Equation (j3.6p for the tensor perturbations are given in [l^. Equation 
(j3.7p for the vector perturbations seem to have been ignored so far. Equation (j3.8p for the scalar perturbations can 
be found in [T7| . 



IV. THE ACTION 



The expansion of the Einstein part (7 = 0) of the action (11.11) at quadratic order in the perturbations on a 
Friedmann-Lemaitre background was first obtained in newtonian gauge in [2ll |. The expansion of the Weyl part is 
easy (see Appendix R)) . The result is (all spatial indices being raised with 6''^): 



with, see e.g. [l9|, [22 



c(T) _ 



K S', 



(s) 



I jd'xa^ [- 



12-H *„ + 25,$„ (2a'^'n + a*$„) - 2{n' + 2%^) 



+ « (Xn - a.Xn - V^^ xl - 6¥>' Xn - 2(^' - 20^ ^>,, Xn)] 



(4.1) 



(4.2) 



Extremisation of S**-"^-* with respect to the tensor perturbations hij readily yields the equations of motion p.6p . 
Similarly the extremisation of S^^^ with respect to the vector perturbations yields the equations of motion p.7p . 

As far as we are aware the fact that the extremisation of S*'^' with respect to the scalar perturbations $„, Xn and 
also yields back the equations of motion p.8p does not appear in the literature, even in the case of standard 
inflation when 7 = 0. This is however the case, as we show it in some detail in Appendix [Bl 



^ The field equations 11.21 1 have ten components. The remaining three (which are the part of the spatial equations proportional to 5;^) 
can be explicitly shown to be redundant. 
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This shows that one can completely fix the coordinate system from start (instead of keeping the ten metric pertur- 
bations plus Sfjj) , obtain the action at quadratic order in terms of seven perturbations only, and still recover the seven 
equations of motion after extremisation, at least when working in the newtonian gauge. Doing so, we do not loose 
any algebraic constraints (or lower derivative equations if 7 ^ 0). There is therefore no need at linear level to keep 
the coordinate system unspecified (that is, keep the lapse and shift as free Lagrange multipliers) as is necessary in any 
hamiltonian formulation of the full theory, and as is usually done in the theory of linear cosmological perturbations, 
see e.g. or [22|. 

V. EVOLUTION 

We shall work in Fourier space, that is, we expand the seven perturbations hij {rj, x'^), ^i{r], x^), <&n(??, *n(?7, x^) 
and XniVi^'')^ collectively denoted by f{r],x''), as 

/(77, x'^) ^ Jd^k 4,(77) '^■^ with r -(r;) = 4,(77) 

and, to simplify notations, we shall omit the index k on the Fourier component fj:{ri). We recall how the Fourier 
components evolve in standard inflation when 7 = in Appendix [U] 

When the Weyl term is present, that is when 7 7^ 0, the structure of the equations of motion (|3.6p p.7p p.Sp 
changes drastically. 

A. Vector perturbations 

Let us start with Eq. (|3.7p for the two vector perturbations ^i. It is no longer an algebraic constraint as in standard 
inflation when 7 = 0. During inflation when spacetime can be harmlessly approximated by a de Sitter space with 
a = l/i—Hr]), Eq. (|3.7p reads, in Fourier space 

with z — —krj. As for ^>i, it stands now for the Fourier component ^E*^ ki''l)- Since (j5.ip is a second order differential 

equation, there are (for each k) two vector degrees of freedom, one for each polarisation i. The two independent 
solutions of (|5.ip for each degree of freedom, that is its two modes, can be given in terms of Bessel (or Hankel) 

functions of index 1^ = y/l/A H^J, cx z^/^ hI''\z) {a = 1, 2). 

Now, 7 must be positive otherwise 4'^ would behave as a tachyon on flat spacetime, see Eq. p.7p . We also have 
that 

^Planck ^ 7 < 

if inflation occurs at GUT scale and if the Weyl-correction is due to a low-energy approximation of some quantum 
gravity theory. 

The evolution of the two modes of each degree of freedom can thus be easily deduced from (jS.ip : at the beginning 
of inflation, when j z^ ^ 1, that is when ^(k/aY ^ 1 (but with ^pianck (^/'^)'^ <C 1 to remain below the 
transplanckian regime), the two modes oscillate with a constant amplitude as in M4. When the "effective mass" term 
1/(7 iJ^z^) becomes dominant the degree of freedom does not tend to a constant (as all degrees of freedom do in 
standard inflation, see Fig. 2] in Appendix ICl) : on the contrary both modes oscillate more and more rapidly as z — >■ 0, 
albeit with a decreasing amplitude: 

^i<^e^" for 7 7/2^2^1; cx Vie=^*^"^/(v^^) for -^H^z"^-^!. (5.2) 

See Fig. m 

Here it may be worthwhile to note that if one considers the case ^ ^ 1, the asymptotic behaviour of the vector 
modes become oc z^^^^"^, where v ~ 1/2 — 1/(7 i?^) as z — ?> 0. 

In the general case now, that is for generic 0(77), equation p. 71) reads 

+ (^e + = (5.3) 
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FIG. 1: Evolution in cosmic time t — j^adrj of a Weyl vector mode 5R[^'j on a de Sitter background (the same behaviour 
is observed in power-law inflation, acct'', p> 1). Values of the parameters: 7 = 1/(25 i?^), k = 2H. 



whose asymptotic zero-mode solutions are given in the WKB approximation by 

cx^e±**/v^ (5.4) 

hence showing the generality of the result (j5.2p . 

In the newtonian gauge identifies to the metric perturbation Bi. In any other gauge we have = Bi — E'^, see 
Eq. (j3.3p . where either Bi or Ei are chosen at will. 



B. Tensor perturbations 

Let us now turn to the Eq. (j3.6p for the two tensor perturbations hij. It is a fourth order differential equation 
which hence describes no longer two, as in standard inflation, but four degrees of freedom (just like in flat spacetime, 
see HJ). During inflation when spacetime can be described by a de Sitter space with a = l/(— r/), the exact 
solution of Eq. (13. 6p is known [l^. It can be written, in Fourier space, as 



Hence the two "Einstein" degrees of freedom h^^' cx e^*^ (1 =F iz), which are the same as in standard inflation, first 
oscillate, decreasing as z, and then tend to constants when z <^ 1, see Fig. 2] of Appendix [C] As for the two Weyl 
degrees of freedom they behave like the two vector degrees of freedom and always oscillate, see Eq. (|5.ip and 

Fig. [T] Hence the amplitude of h'^^^ = /l^J^-'/a, first decreases as z, and then as z'^/^ as 2; — > 0. Therefore the full 
Fourier component of the metric perturbation hij first oscillates and, as infiation progresses, the standard Einstein 
modes eventually dominate, tending to a constant, see Fig. [2l 

In power-law inflation, a{t) oc with p > 1, where t is cosmic time, the equation of motion p.6p does not split into 
two second order differential equations, one for the Einstein degree of freedom the other for the Weyl ghost, as when 
the background is de Sitter spacetime. Its solutions however behave similarly, see Fig. [51 

More generally we find that the two Einstein zero-modes, solutions of h'i'j + 27i hij = are approximate solutions 

of the full zero-mode equation of motion, h'-j + 2% h'ij + 70"^ = 0, if j H"^ <C 1 (with H — H/a). As for two 
Weyl zero-modes they can be found using the WKB approximation so that, all in all, the four independent tensorial 



(E) 



-(E) 



(W) 



7^2 22 



(5.5) 
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FIG. 2: Evolution in cosmic time of a Fourier component of the tensor metric perturbation hij on a de Sitter background (the 
same behaviour is observed in power-law inflation). Values of the parameters: 7 = 1/(25 f/^), k — 2H. 



zero-modes behave as 



dt 
7^ 



hence confirming that as infiation progresses a generic linear combination of the four modes will tend to a constant. 

C. Scalar perturbations: a master equation and its solutions 

The analysis of the three Eq. (|3.8I) for the scalar perturbations is slightly more involved. However it is easy to 
extract from them a master equation for W , which reads, in Fourier space: 



7 ( VF*") - ) + C2 ly + Ci -t- Co = , 



H 



(5.7) 



where W = ^'n — <&n7 where a dot means derivation with respect to cosmic time t — J^a drj, where H = d/a, and with 



Cn = 2H- 



' \^ ~ 




H 

^ +7 
H 
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Hii 

— ^ -+ 


— + 


H 





fc2 IH_ 
a? \ H 



AH 



(5.8) 



H 



- -H + 2H^ + ^ . 



(When 7 = 0, we have that W = — 2<i>n (see Eq. (j3.8l) ) and Eq. (|5.7p reduces as it must to the equation for $„ in 
standard infiation given in Appendix [CJ Eq. (|C4p .') 

Equation (j5.7l) is fourth order and therefore describes two degrees of freedom, and not only one as in standard 
inflation (see Appendix ICj) . Ideally one should try to decompose (15. 7p into two second order differential equations 
for an "Einstein" and a "Weyl" degree of freedom (as was done in 16] for the tensorial perturbations on a de Sitter 
background, see above). We leave this to further work and content ourselves here with an analysis of the solutions of 

. . . 

To have a grip of their behaviour we assume power-law infiation: a{t) oc with p > 1. 

Proceeding along the line that we followed to analyse the tensor modes we find that the standard infiationary 
zero-modes which solve Eq. (j5.7l) when 7 = approximately solve the full equation if ^ H"^ <C 1. As for the other two 
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FIG. 3: Evolution in cosmic time of a Fourier component of the scalar mode W in power-law inflation. Values of the parameters: 
a = ao (t/toY, Ho = p/to with p = 4, and 7 = 1/(25 Hq), k = 2Ho ao; initial conditions at t ^ to: W ^ 0, W ^ 1/to, W ^ 0, 
W'^'' = l/tl 



modes they are found using the WKB approximation, so that the four independent zero modes of Eq. (j5.7p have the 
following late time behaviour 

W(x{l, i5e"/v^, t^e-'*/^}. (5.9) 

(This is confirmed by an analysis of the leading behaviour of the solutions at the irregular singular point at infinity, 
as well as the exact zero mode solutions which can be written in terms of Bessel and hypergeometric functions.) 

These behaviours are in striking contrast to those of the tensor modes which are dominated by the constant, 
Einstein-mode. Here both Einstein modes are subdominant. 

The evolution of a typical Fourier component of W is given in Fig. [31 As one can see, not only does the Fourier 
component W never "freeze out" but its amplitude increases as inflation proceeds instead of tending to a constant as 
in standard inflation when 7 = (see Appendix [CJ. 

Equation (j5.7p and the behaviour (j5.9p of its modes can be seen as the main result of this paper since, knowing W, 
we can predict the behaviour of all cosmological perturbations of the model. 

D. Evolution of the scalar perturbations in the newtonian gauge 

Once W is known, $„ and Xn follow from Eq. 



-2$„ = W + -f [W + HW+ —W 



KipXn = W + HW + -f 



k 

W^(3) + w {H - H^) + ^{W - HW) 



(5.10) 



It follows from (j5.9p and (|5.10p that ^'^ increases as W and that $„ and Xn behave as W^-^ e±'*/V^. (The leading 
term in $„ grows a priori like W but cancels out; as for tpxn oc Xn/t it should also a priori grow like W but the two 
first leading orders cancel out.) 

Since <J>n = C, ^'n = A and Xn — in the newtonian gauge, we therefore reach the conclusion that in that gauge 
all cosmological perturbations blow up in single field inflation with a Weyl term. 

This result does not mean however that they blow up in all coordinate systems, as we see now. 



E. Evolution of the scalar perturbations in the comoving slicing 



All linear combinations of $n, ^'n and Xn are gauge invariant and can be expressed in terms of W only. Moreover, 
using Eq. p.3p . they give the perturbations {A, B, C, E) of the metric and the perturbation 6(j) of the scalar field in 
any gauge. Therefore gauge invariant quantities can be built which identify to various perturbations in a given gauge. 
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An example a such a gauge invariant perturbation is the curvature perturbation [26 



7^e^<i>„-^=C-^ (5.11) 
ip' If' 



(introduced in standard inflation in various guise, see [l9|, |26|, |27| and Appendix ICj) . In the "comoving slicing" gauge, 
that is, in the coordinate system where 

(50 = 0, (5.12) 

Tic (denoted by —C, in [1^) identifies with the metric perturbation C. 

From its definition (|5.1ip one expects a priori that TZc should grow like 3>n, that is as t^l'^~^ g***/^ or perhaps at 
a slower rate if the first leading terms cancel out, but it happens that the cancellation is so drastic as to make 7?,c not 
to grow (despite the fact that p can be as big as one wishes). 

Indeed, 7?.c, using (I5.10p . is a function of W up to W^"^^ . If we compute Tic and use the master equation (j5.7p to 
eliminate VF*^'*-' we find that 



2a2 1 ii 



3 H / r\ a? 



(5.13) 



(When 7 = we recover the well-known result of standard inflation, see Appendix [Cl) 

The analytic analysis of TZc, as well as numerical plots show without any ambiguity that it goes to zero as 
^-3p/2 gilt/ s/i _ must therefore conclude that indeed all leading terms up to order t~3p/2 ^±it/ ^ cancel out 
when looking carefully at the asymptotic behaviour of 7?,c from its deflnition. More precisely, knowing the asymptotic 
behaviour of the four independent solutions for VF, see (|5.9p we have that 



Consider now the following gauge invariant perturbation: 



which identifies to the metric perturbation A in comoving slicing. 

Again one expects a priori from its definition that A^, will grow, like ^'n that is like W if the leading orders do not 
cancel out. But, again, this not the case. Indeed A^ depends, using (15.101) . on W and its derivatives up to the fourth. 
Using the master equation (|5.7p to eliminate W'^^^ we have that 



Ac 



2a? E 
We also have that 



-i{W - HW) + W I 1 + ^ 



(5.15) 



nc-HAc = ^W. (5.16) 

Knowing the asymptotic behaviour of VF, see (|5.9p . we see that the growing modes solve the equation 7 (W^ — 
H VF) -t- = 0. Therefore the W mode which gives the asymptotic behaviour of A^ is the subdominant mode W —\. 
Hence A^ decays to zero as without oscillating (contrarily to T^c), as numerical plots confirm. 

We therefore have shown that the metric perturbations A and C, do not grow when evaluating them in the comoving 
slicing. This is in striking contrast to their behaviour in the newtonian gauge, where they blow up, as we have seen 
in the previous section. 

A last check has to be done though, since in the comoving slicing the metric perturbation {B — E') is not zero. To 
find it one uses the expression for in terms of Scf): Xn = 54> + f' {B — E'), see Eq. (|3.4p . In the comoving slicing 
S(l> = 0. Therefore B — E' = Xn/y^' = Xn/iatp). We know from the previous section the asymptotic behaviour of Xn^ 
Xn oc e±^*/v^. Therefore {B - E') also decays, as t-P/'^ e±**/V7. 
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VI. CONCLUSIONS 

In this paper, we have studied the role of the Weyl term on the evolution of linear cosmological perturbations when 
the Friedmann-Lemaitre background is that of single field inflation. 

We found that the two Weyl tensor degrees of freedom are tamed by the Einstein gravitational waves and thus do 
not spoil too much the evolution of the tensorial cosmological perturbations as given by the standard inflationary 
scenario, whether the background is approximated by a de Sitter spacetime as in [l6l | or in power-law inflation, see 
Fig. [2] Vector modes on the other hand, which are absent in standard inflation, do propagate; these two pure- Weyl 
vector degrees of freedom never "freeze out" , but their amplitude decreases as inflation proceeds, see Fig. [T] Finally, 
the evolution of the scalar modes is drastically modified by the presence of the Weyl term: instead of one there are 
now two scalar degrees of freedom which, when working in the newtonian gauge, and contrarily to what happens in 
standard inflation, not only do not freeze out but their amplitude increases during inflation, see Fig. |3l However there 
exists at least one coordinate system (the comoving slicing) where none of the perturbations grows. 

We cannot therefore claim at this stage whether the five Weyl degrees of freedom, which are ghosts in Minkowski 
spacetime, screw up or not the evolution of linear cosmological perturbations in inflation, since their asymptotic 
behaviour depends crucially on the coordinate system used. 

To complete our study, and arrive at a more definite conclusion, an hamiltonian analysis of the action (|4.1| - |4.2|) 
should be performed to isolate the Weyl degrees of freedom from the Einstein's ones. It is however clear from the form 
of the action that the vector perturbations are two ghosts, since the sign of their kinetic term is positive; it is also 
clear that their quantisation will impose a normalisation of their Fourier modes in fc^"^/^ because of the presence of 
the extra spatial derivatives in the kinetic term in the action S*^^^ . As for the tensor perturbations they were analysed 
in [i^ when the background is approximated by a de Sitter spacetime: the two Weyl degrees of freedom /i^J^^ are 
ghosts, and the normalisation of the Einstein modes is modified by their presence. It remains however to generalize 
this analysis to the case when the background is no longer de Sitter spacetime. Finally, the hamiltonian analysis of 
the action 5^^-' (|4.1H4.2p for the scalar perturbations, in order to isolate the ghost degree of freedom, is more tricky 
and is left to further work, see (23 |. 

To complete our study the question of what happens at the end of inflation should also be addressed. As a first 
step this transitory period could be modelled by a sudden transition from the inflationary stage with a (x with 
p > 1 to the radiation era a cx t^^^. The junction conditions which give the perturbations after the transition in 
terms of their behaviour during inflation are well-known in Einstein's theory, see [2^ . When the Weyl-term is present 
they have to be analysed anew since the equations of motion become fourth order. One expects however that the 
two tensorial Weyl ghosts will not change too much the standard picture since they are subdominant compared to 
Einstein's gravitational waves. The matching of the two ghost vector degrees of freedom, although decaying during 
inflation, may on the other hand be more tricky as they oscillate at very high frequency at the end of inflation. Finally 
a proper matching of the two scalar degrees of freedom requires first the hamiltonian analysis referred to above. 

Last but not least a complete study of the role of Weyl's ghosts in inflation requires an analysis of observables, such 
as the CMB temperature fluctuations, which may be affected by their presence. 

In any case we have already seen in this paper that the addition of the Weyl term to the action of Einstein's gravity 
coupled to a scalar field modifies drastically the evolution of perturbations in inflationary cosmological models and 
we gave some of the necessary tools to assess their influence in observational cosmology. 
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Appendix A: Perturbed quantities 

1. Perturbed Weyl and Bach tensors (See definitions in footnote [T|) 

• Linearised Weyl tensor on M4: 

- Scalar part (with W ='ifn- $11): 



C. 



0, 



(s) 

(S) 
Oijk 

^ifki = 7:{^ik dji - 6a djk - Sjk da + 5ji dik) W - \{5ik Sji - 5a 5jk) AM^ . 



(Al) 



Vector part: 



^Oijk 



^ijkl 

— Tensor part: 



-(a.*, + a,*.o, 

^d,{d,^k - dk^j) - i(<5,, A^k - s,k Ai-,) , 
1 



ClJki = i[S^k [dj^i + di^j) - 6a [dj^k + dk^j) - S^k (d.^i + m^) + S,i (d.^k + ^k^^ 



(T) 
OiOj 

(T) 
Oijk 



(T) 
ijkl 



^-^djhik + -^dkhij , 

-{-dikhji + dahjk + djkha - djihik) + ]r{Sik Ohji - 5a ^hjk - 5jk Uha + S^i Uhik) 



(A2) 



(A3) 



Since the perturbed Friedmann-Lemaitre metric is conformal to the perturbed Minkowski metric, and as can be 
checked exphcitly, the components of C^^pa are the same for both metrics. 

• Linearised Bach tensor on M4: 
At hnear order around M4 the Bach tensor reduces to Bf^^ ~ d'^'^Cf^pi^a- Its components are: 



B, 



(S) 

00 



1 



-AAW^, B, 



(S) 



1 



^a^Aiy', B^^ 



(S) 



^a., (w" - ^AW] i% iW" AW) , 



(V) 



0, B, 



(V) _ 



--AD^i, B 



(V) ^ 



(T) 
00 



0, i3^P=0, BiJ^^-^ODh,,. 



(A4) 



The hnearised components of _B^^ on a Friedmann-Lemaitre metric are the same, up to a factor 1 /a^ 



Weyl action at quadratic order 



Using the expressions for the hnearised Weyl tensor on M4 given above and because of the conformal invariance of 
the Weyl action we have 



^ I d X (-\/ ~QCpi/p(j ^ )pcrturbcd FL — ^ j d X (\/ ^Cpi^pcr ^ )pcrturbed Minkowski 



= I d^X 



{AW)^ + 7(9,;*j d'i/^ - A^i A¥) + li'h.j h'^ - 2dkh,j d^h'^ + A/i^ Ah'^) 
4 8 



(A5) 
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Appendix B: The equations of motion from the action for gauge invariant perturbations 

NB: in this appendix we suppress the index n which ornate $ and 4* and x i^i the main text. 
Consider the action S^^'^ in terms of the three perturbations $, ^' and x given in (j4.ip (|4.2p . 

Extremisation of S''^^ with respect to 'I' gives the (OO)-component of the perturbed equations of motion, already 
obtained in p.Sp . that is: 

m $' - 2A$ - (ip'x' + 2a^V'^ + x) = ^ AA(* - $) . (Bl) 

Extremisation of S'^^^ with respect to x gives the perturbed Klein-Gordon equation 

X" + 2H x' - Ax + a' V^^ x-v' - 3$') + 20^ 1/^ ^ ^ (B2) 
Finally, extremisation of S*'^-* with respect to $ gives 

A f* + $- -^AC*-*)) = 4r(a'i^)' with + (B3) 

Let us now perform the following manipulations: 

1. Replace in (IBip x ^-nd x' by their expressions in function of F given in (jBSj; (|B1|) then depends on 5', $ and F. 

2. Replace ^ by its expression in terms of F and $ using (IB3p : (IBip then depends on $ and only. 

3. Extract from it the expression of $" and compute its time derivative 

Consider now Eq. (|B2p : after the above replacements of x ^-nd it depends on F, $ and its time derivatives up 
to the third. Replace and by the expressions obtained previously and find (we used Mathematica) that it can 
be written as 

AG = with G= ^^^^ + ($'--H*) + -^A(«'-$)'. (B4) 
2 3a^ 

Therefore the Klein-Gordon equation (jB2l) is equivalent to the (Oi)-scalar component of the equations of motion p.8p . 

As for (jB3l) . since G = and hence F = — g^A('I' — <&)' with F given in (|B3p . it becomes the (ij)-scalar components 
(i ^ j) of the equations of motion (13. 8p . 

This (simple) derivation must be contrasted to those found in the literature. In [l^ and [l^l for example, KS'i'^ is 
computed without specifying the coordinate system, that is for the full metric 

ds^ = a{-qf {-(1 + 2A) drf + 2d^B dx' drj + [S^-j (1 + 2G) + 2d^.jE] dx' dx^} 

and the action (in the case 7 = 0) hence gets an extra-term, proportional to {B — E') F, with F = '^'^ ^ + ($' — "H ^'). 

Extremisation with respect to {B — E') gives F = that is, the (Oi)-component of the equations of motion given 
in (13.81) (for 7 = 0). Extremisation with respect to $, Eq. (jB3l) . then gives the last equations of motion p.8p . (And 
the perturbed Klein-Gordon equation is ignored.) 



Appendix C: The evolution of cosmological perturbations during standard inflation: recap 

In the case of standard inflation (7 = 0), Eq. p.7p is a constraint and vector modes are absent: 

^,=0. (CI) 



Equation p.6p describes the two degrees of freedom of gravitational waves that freely propagate on the Friedmann- 
Lemaitre background, and can be rewritten as (23j : 

a" _ _ _ - 

n'lj —^ij - Afiij = with fiij = ahij . (C2) 

As for Eq. p.Sp . it consists in two constraints: 

= and 1^=n^r.-K (C3) 
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FIG. 4: Typical evolution of a cosmological perturbation mode "I>n during inflation in Einstein's theory. 



which, inserted into the first of Eq. p.Sp , yield a master equation which can be written into various equivalent forms 
[EllflllM, e.g.: 



2 n 



$' + 2 -H' - -H 



that is, 



(C4) 



a H 
-u — Au = with u = — $n and 9 = ; 



Using the constraint (jC3p this master equation also reads 



V V — Av = with z : 

z 



a (fi 
IT 



and V = a \ -j^^n 



Xn 



2a 

K if' 



a tf' 



(C5) 



(C6) 



The evolution of the two degrees of freedom fiij and the evolution of u (or v) are similar in the inflationary stage 
when the "mass-terms" a" /a, 9" /9 or z" / z can be neglected: in Fourier space the two independent modes of each 
degree of freedom oscillate as in flat spacetime. When inflation progresses and the mass-terms come to dominate 
only the dominant modes of hij and ^>n "survive" and become almost constant, see Fig. |4l More precisely the two 
zero-modes which solve (|C5p behave as 9 and 9 J^drj/9. In the case of power-law inflation, a oc this translates as 



$„ cx {1 



-(i+p) 



}• 



(C7) 



After the end of inflation when the scale factor increases as t^/^ say, the mass-terms become subdominant again and 
the modes again oscillate. 

Of course other gauge invariant variables can be introduced which, thanks to the constraints (IC3|) . can all be 
expressed in terms of $n- An example is the curvature perturbation [26|: 



ifi' Z 



2n 



K if 



/2 n 



(C8) 



(denoted by — C in ilSi])- It is a useful quantity because its time derivative is given by (using the master equation (|C4 
to eliminate $"): 



(C9) 



Therefore the amplification of the $„ modes can easily be obtained from the fact that TZc is almost constant as long 
as the mass-term dominates. 
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Another interesting gauge invariant perturbation is 



1 / oXi 



) 




2a 



II 



(CIO) 



n 



a\ If' 



a \_Kip' 
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which, using (jC4l) . can be shown to be simply related to TZc by: HAc = Tlc. 

All gauge invariant variables can now be related to the perturbations of the metric and the scalar field using 
$„ = C + 'H{B — E') (Eq. p.3p in the main text). Thus, in the newtonian gauge, where B — E — the gauge 
invariant perturbations $„, and Xn identify respectively to C, A and Scj). 

As for TZc and A^ they are related to the metric and scalar field perturbations as 



In the "comoving slicing", that is, in the coordinate system where Scf) — 0, TZc and Ac identify to the metric pertur- 
bations A and C. 
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